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I. INTRODUCTION 

The holographic principle , especially with its first realization in string theory-the 

AdS/CFT correspondence, offers us very intriguing and powerful tools to deal with the strongly 
coupled quantum systems from the dual viewpoint Q-S]. The more general framework of the cor- 
respondence, which is called the gauge/gravity duality, has been extensively applied to the study 
of the QCD, Quark Gluon Plasma, hydrodynamics etc., for an incomplete list, see 6h28||. In the 
framework of the gauge/gravity duality, the features of strongly coupled quantum field theory on 
the conformally flat boundary can be fully captured by its dual weakly coupled classical gravita- 
tional or string theory in the curved bulk spacetime. Even though the gauge/gravity duality is 
widely believed to be held for arbitrary spacetime backgrounds, so far there are only a few explicit 
examples, in which the best known one is that the strongly coupled A/" = 4 supersymmetric Yang- 
Mills theory in four dimensional flat spacetime is equivalent to the classical (weakly coupled) limit 



3 



of the type IIB superstring theory (supergravity) in AdSs x spacetime. For most other cases 
one still requires the bulk to be asymptotically AdS spacetime whereas the boundary field theory is 
conformally invariant and relativistic. However, besides numerous strongly coupled systems in high 
energy physics described by the relativistic quantum field theory, there also exist large classes of 
strongly coupled phenomena described by non-relativistic field theory in various condensed matter 
systems, especially near the (quantum) critical points. Therefore, it is very interesting and impor- 
tant to extend the gauge/gravity duality into a non-relativistic version in order to understand the 
strongly coupled phenomena in the laboratory conditions. 

Much progress has been made towards this direction in the past few years. One class of work 
focused on the study of field theories with the Schrodinger symmetry, motivated by the study of 
fermions at unitarity, see 2^, 3^. Another class of work tried to utilize the dual gravitational 
theories to study the condensed matter systems near quantum phase transitions that contain the 
Lifshitz fixed points 
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such as the strongly correlated electron systems. The particular 



property of the Lifshitz symmetry is that it consists of the anisotropic scaling 



Xx and t — )• X^t, 



(1) 



where z is called the dynamical critical exponent. When z = 1, the above transformation is the 
usual relativistic scaling. From the perspective of the gauge/gravity duality, the essential point is to 
construct bulk gravitational solutions by adding some appropriate sources to realize the boundary 
non-relativistic quantum field theories with the Lifshitz symmetry. The first attempt was done in 
3ll |. in which a four dimensional asymptotic Lifshitz spacetime at zero temperature was obtained 
in the AdS Einstein gravity together with 1- and 2- form gauge fields. The bulk solution can be 
viewed as a toy model to provide us some useful descriptions for certain magnetic materials and 
liquid crystals. Subsequently, many asymptotic Lifshitz black hole solutions have been found and 



analyzed, see for example 
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36|, 



42 



47|. With the help of these solutions, important properties of 



the dual strongly coupled non-relativisitc quantum field theories, such as the transport coefficients, 
n-point correlation functions, renormalized stress tenosr and higher order corrections 48145 can 



be studied by performing the calculations on the side of the Lifshitz black holes/branes. 

The asymptotic Lifshitz solutions can be obtained from different types of theories, the one 
received much attention is the Einstein-Maxwell-dilaton (EMD) theory, which may be consistently 
truncated from higher dimensional supergravity theories. Recently, a class of analytic Lifshitz 
black hole/brane solutions have been solved in the EMD theory by adding multiple independent 
C/(l) gauge fields 64]. These kinds of charged Lifshitz black hole configurations can provide 
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potential interesting applications to condensed matter systems such as fluids, non-Fermi liquids 
and conductors that contain the Lifshitz fixed points. Some holographic aspects in these spacetime 
backgrounds have been brought out, such as the instabilities of dual superfiuid by adding probe 



charged scalar field in the bulk 



I52I ti 



5514571]. 



to further study 



53( 1 ^. For other related works, see for example 
The purpose of this paper is to utilize these charged Lifshitz black branes 
certain interesting phenomena of the dual strongly coupled non-relativistic quantum field theory 
with the Lifshitz fixed points on the boundary. Based on the dictionary of the gauge/gravity 
duality, we know that the multiple U{1) gauge fields in the bulk will source multiple electric 
currents in the boundary field theory. As a theoretical model, there is no constraints on the 
number of independent electric currents even though their physical interpretations are not yet 
very clear. What we focus in this paper is to investigate the transport coefficients of the dual 
non-relativistic field theory, which includes the electric conductivity cr, the thermal conductivity 
K, the thermoelectric conductivity a and the shear viscosity i]. To reach this goal, we consider the 
linearized gravitational and gauge fields perturbations (the scalar channel and the shear channel) 
in the bulk EMD theory. In particular, the bulk Lifshitz black hole can be viewed as the non- 
relativistic counterpart of the Reissner-Nordstrom-AdS black hole when N = 2. Focusing on this 
case we calculate the conductivities of the dual non-relativistic field theories numerically, which 
are expected to capture the universal behavior of a class of conductors near the Lifshitz fixed 
points. Speicifically, after deriving the renormalized second order on-shell effective action, we 
work out the numerical results of conductivities, including the electric, thermoelectric and thermal 
conductivities. In particular, we work in d = 3 and d = 4 (d is the dimension of the boundary 
field theory) for 1 < 2; < 2. We find some new frequency dependent power law features of the AC 
conductivities in the large frequency regime for 1 < z < 2. The possible relations between these 
novel features and the 'symmetric hopping model' in condensed matter physics are discussed in 
the context. In addition, an other interesting problem is to see whether these additional bulk U{1) 
gauge fields and dilaton will affect the famous Kovtun-Starinents-Son (KSS) bound derived in the 



Einstein gravity 
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ll[ |. By solving the equation of motion of the transverse graviton at the low 
frequency limit and applying the linear response theory, we show that this bound is not violated 
although the additional gauge fields and dilaton do respectivley contribute to the shear viscosity 
as well as the entropy density of boundary charged fluids. 



^ A generalization of theses solutions with additional hyperscaling violation factor was obtained in [Hil, in which 
their dual nonrelativistic field theories was briefly analyzed as well. 
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The outline of the paper is as follows. In Section [iTl we give a brief review of the Lifshitz black 
hole/brane backgrounds that we will use in this paper. In Section Hill we obtain the renormalized 
second order on-shell action of the perturbations and compute the electric, thermal and thermo- 
electric conductivities of the boundary non-relativistic field theory, in the N = 2 case. We calculate 
the shear viscosity of the boundary fluid both for N = I and generic N cases by solving the equa- 
tion of motion of the transverse graviton in Section llVi Conclusions and discussions are drawn in 
Section |Vl Besides, we list some detailed calculations for deriving the perturbation equations and 
the second order on-shell actions in Appendix A. 



II. THE CONFIGURATION OF LIFSHITZ BLACK HOLES/BRANES 



Let us consider the {d + 1) dimensional EMD theory 



dV{(l)) l^d'Ja 2 



I = J d^^'xV^ - - -{d^r - -J{4>)A' - y(c/.)J , (2) 

where F^i, = d^Ay — d^A^ is the U{1) gauge field strength, (p is the dialton field, 7((^) is the 
coupling between the gauge field and the dilaton, J{(t>) is the source term and V{(j)) is the potential 
term. When we add > 1 number of independent U{1) gauge fields, the above variables can be 
accordingly changed as "f{(j)) — )■ Yla=i la{4') and F^jy — )■ Fa^u, and the equation of motions are 

N 

v,{U4>)Fr) = ja:, 

1 1^/ 1 

a=l ^ ^ ^ 

+i (^d^(Pd,(t> - \gAd<t>? - g^..v{(t>)^ . (3) 

The Lifshitz black holes can be obtained from the following ansatz 

ds^ = -^(r)e-^^''Ut^ + — + b'^{r)dx'dxi, (4) 



together with 



Aa = Aat{r)dt 



(r), J(</.) = 0, y((/>) = 2A and 7a = e^'^'^'. (5) 
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For = 1 case, the solution is 



H(r) = 1 - A'^t = l-'^j2{d + z-l){z-l)iiy ^^(^r'^+^-s, g'^ = /irV2{--i){'^-i), 

/ 2(d-l) (^ + rf-2)(^ + d-l) 

^ = -V^^' ^ = ' 

where / is the curvature radius of the Lifshitz spacetime, ^ is the scalar field amplitude, m is related 
to the mass of the black hole and ""' is the derivative with respect to r. The Hawking temperature 
and the Bekenstein-Hawking entropy are respectively 



{z + d-l)rl Vd^frhY-' 
and Vd~i = f d'^~^x is the spatial volume of the boundary. 



(7) 



For generic A^, the black hole solution is 



52| 



2z j2 2 



f,(r) = k( ?— + 1 - mr-('^+^-i) + V - r-'^{d+z-2) 



d-l 



z-l 



= p,/.- V ^^^ir2-^2^Ja^2_ . , A^ - 1) 

/2(d- 1) , /2(z- 1) , d-2 /2((i- 1) , ^ ^, 

where pa are related to the charges of the black hole, while k is the factor indicating the topology 
of the horizon. For A; = 0, the horizon is flat; for k = —1, the horizon is hyperbolic and the horizon 
is spherical for k = 1. In the following, we shall take the spatial flat case, namely, the Lifshitz black 
brane with A; = 0. When N > 2, the black brane will contain multiple horizons in the presence of 
electromagnetic fields, let us define the outer event horizon to be located at r = r/j, i.e. f{rh) = 0. 
Then the temperature of the black brane is 

T 1 rM'+'fV ^ 1 /r;,y+i /2(d + z-2) m{d + z-3) \ 
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where /(r) = 1 — j77,7--('^+2-i) -|- J2f=2 2(d-i)(d+z-3) r~^^'^~^^~'^^ . The entropy Sbh and entropy 



density s of the dual CFT are 



III. THE CONDUCTIVITIES 

In this section, we will compute the conductivities of the non-relativistic quantum field theory 
dual to the Lifshitz black brane. The electric conductivity a can be calculated by just turning on the 
bulk gauge field fluctuations 5Ax{t,r) = ax{r)e~^^^ . However, if we want to consider the thermal 
conductivity R and the thermoelectric conductivity a , we need to consider the back reaction of 
the gauge fields to the metric, namely, we need to meanwhile turn on 5gtx{t-,r) = htx{r)e~'^^^ . For 
the EMD theory (when taking J = 0) in eq. ([2]) , we can obtain the linearized Einstein and Maxwell 
equations as (see Appendix A for details) 

26' ^ 

Kx - -^^te + X] 7a(0)^ataax = 0, (II) 
a=l 

j^dj^\ _(2_y}Hx_Kx\ 

where 6(r),^(r) and x(r) are factors in eq.([l]). Note that eq. ()lip is the first order differential 
equation for htx which can be integrated out as 

1 ^ 

htx = -h{rf j ^Y/1,{(t))^^^aaxdr, (13) 
and eq. p2p can be written into the following equation as 



^ax + I ^ + - - — + 



) o^'ax + ^e>^aax = ^ {^lc{<P)a,xA'}l ^;,e^ (14) 



b e 2 -fa{^) # 

with the help of eq. (|ll|) . 

When = 1, the background Lifshitz black brane eq.® is neutral as the Schwarzschild AdS 
black brane, the electric conductivity has been studied by adding a probe U{1) gauge field in the 
bulk in 0- 

In the following, we will focus on the N = 2 situation, in which 

e-^=(y)''"'' e(r) = ^/(r), = ^rV^C'^-DC-D, b{r) = j, 

/(r) = 1 - mr-^'^+^-i) + , f^^ ' , -r-'('i+^-') . (15) 
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Recall that for the N = 2 case, the background gauge field An is divergent at the spatial infinity, 
it only supports the asymptotic Lifshitz geometry instead of contributing to the free charge of 
the background electromagnetic field 



52l |. On the contrary, the gauge field A2t plays the role of 
the free electromagnetic field. Besides, our numeric results show that the asymptotic expansion 
of aix is also divergent at the spatial infinity. Thus only the fluctuation of A2, namely, a2x is 
the genuine electromagnetic perturbation, which will contribute to the conductivities of the dual 
field theory on the boundary. Consequently, to study the conductivities, we only need to turn on 
the perturbations a2x and htx, while turning off the perturbation aix- Then after substituting the 
above black brane solution eq. (ll5[) into the fluctuation equations (llip and (I14p . we obtain 



Kx htx+ P2r 

r 



z—d^ 



«2x + ( y + ; J a2x + 



'j^^{2-2d~2z)i2z 



f 



a2a 



022 



0, 



0. 



(16) 
(17) 



The explicit asymptotic behavior of a2x near the infinite boundary with certain d and z considered 
in this paper can be found in Table HI in which Ci and C2 are expansion coefficients that depend 
on the frequency uj. According to the gauge/gravity duality, Ci represents the source while C2 
represents the vacuum expectation value of the current operator Zx dual to a2x- 





z = 1 


z = 3/2 


Z = 2 


d = 3 


Ci + ^ 

1 r 


Ci + ^2 


'-'1 "t" "T 7^ 


d = A 




n 1 2CxiJ- , C2 


n 1 C\J^ 1 C2 



TABLE I: The expansions of aix with respect to various d and z near infinity. The coefficients C\ and C2 
are functions of the frequency uj. 



In addition, the asymptotic behavior of htx near the infinity boundary is. 



j,{d—z~\) 



+ ••• 



(18) 



where, = Ci/32/(l + d — z) in which C\ is the source term of the expansions in a2x, see TableHl 



Second order on-shell action 



In order to compute the transport coefficients of cj, a and k, we need to know the quadratic 
on-shell actions for these perturbations. The on-shell action for the perturbation a2x and h^x up 



9 



to 2nd order is (we have set / = 1), 



c(2) _ o(2) , M 
"-"on-shell ~ '~'a2x ^ "^htx' 



:(2) 



where 



5(2) 



1 



d-3 



5 



(2) 



x') 



(19) 

(20) 
(21) 



Usually the on-shell action eq. (jl9p is divergent near the asymptotic boundary, the divergence can be 
eliminated through the holographic renormalization approach, i.e. by adding appropriate boundary 



counter terms to the action (see, for example |58l. l59l|). In the configuration of the Lifshitz black 



brane, the counter terms have different forms with respect to different d and z. We will list them 
in the following: 

First of all, we will introduce the counter terms to Sa2i in eq. (l20p . These counter terms are 
classified according to the expansions of a2x in Table HI 

(a) , d = 3,z = 1 and 3/2: 

(2) 

In this case, the on-shell action of ^aj^; is finite at the infinite boundary. There are no counter 

(2) 

terms to Sa2i just like in the usual relativistic holographic superconductors [6C|. 

(b) . (d = 3, z = 2) and (d = 4, z = 1) : 

(2) 

For (d = 3, z = 2) and (d = 4, z = 1), there will be logarithmic divergence for Sa2l on the 
infinite boundary. In this case, the generic expansions of a2x near r — t- 00 now is, 

"2^(^) ~ + (d + 3z-5)r^+3.-5 +g2 [-) 

(2'] 

The divergent term of Sa2i can be obtained from eg. ([20]) as, 



(22) 



I. 



div.a2x 



%iCi2a;2log(r)/iA 



(23) 



where T = 1//3 is the temperature of the boundary field theory and -^p^ is just the volume 
integration J d^x. Therefore, in this case the counter term should be, 



^ct.a2x 



-log(r) / d'^xV^^^mF'?- 



(24) 



where, 7*^ is the determinant of the induced metric while F-j is the induced gauge field strength on 
the asymptotic UV cut off boundary, respectively. It is easy to get that {F-'-)'^ = 2ijJ^{a2xf"r~'^^ j^- 



Therefore, the finite on-shell S'qjI is. 



l^l = Sfl + = / d'^x {c,C2{d + 3z - 5) - ^^^) 



(25) 
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(c). d = 4,z = 3/2 and 2: 

For d = 4, z = 1 and z = 3/2, the general expansions of a2x is, 

«2.(r) ~ + (iT + cJ iT^' . (26) 



2(z - l)z \r 

(2) 

In this case, the divergent term of Saai is , 

The counter term for this divergence now is, 

W = / d'^xy^ l2mF?,?- (28) 

Therefore, from the expansions, we can get the finite on-shell action as, 

Ifl = Sfl + ht.a,^ = I d'x C,C2{2>z - l)^^/^(^. (29) 

Next, we will introduce the counter terms for the on-shell action S\^^ in eq. ipT]) . We can expand 
it near r — )• oo as. 



^ht^ ~ -^div./it, + hnitc.htx 1 (30) 



where. 



Wh.. = j d'x{ht^f{z-2y~^^\ (31) 
/finite.... = / d'x ( ^;j!^;f giMg + lih'H^fmid + z - l)r2"2^) . (32) 

It can be found that when z = 1 the last term in /finite. hta: finite while for z > 1 it will vanish at 
r — )• oo. As usual, we can introduce the Gibbons-Hawking term Igh and a counter term for the 
cosmological constant Ict.cc into the on-shell action to cancel the divergence, where 

Igh = 2 j d'^xy^K, (33) 

Ict.cc = 2 [ d'^xy^id-l). (34) 



in which, K = 7^jyV'^n'^ is the trace of the extrinsic curvature while n'^ is the outward pointing 
unit normal vector on the boundary. Expanding eq. (|33p and eq. (|34p to the quadratic order of the 
perturbations near r — )■ oo, we arrive at, 

igl = I d'x {{hfjfiz - d - ly^^^' + l{h^J)Md + z - ly-'^^ , (35) 

/il = Jd'^x l^ih^y{d-iy~^^' + ^^^^C,p2h^^ + (36) 
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Therefore, the total finite on-sheh action of the perturbation htx can be obtained from eqs. (|3ip . 
(I32D. (I35D and (ESD as, 



4'^ = 4'^ - iSl - iSl (37) 

hta; tltx '-'H Ct.CC \ I 

(-C„,C-=<;^(ft£))V--) (38) 
-f'x 1 - - ^^^CSVr'-^') . (39) 

So, finally the total renormalized quadratic on-shell action for the perturbations 02a; and ht^ is, 
41 = 42 = / + 3z - 5) - (d + 1 - z)}^}^ - ^^^(^1?) V--) , (40) 

for (i = 3, 2; = 1 and 3/2; Or, 

j(2) ^ j(2) .(2) 
total a2i ^ feti 

= /d'^x (w2(d + 3.-5)--^-^-(d + l-z)/.L%«-^^^ 
when (d = 3, 2 = 2) and (d = 4, z = 1); Or, 



/(2) = j(2) , j(2) 
total ci2x ^ -'ft.tj; 



.(42) 



when d = 4, z = 3/2 and 2. 

B. The electric, thermoelectric and thermal conductivities 

As long as we get the quadratic on-shell action for the perturbations, we can derive the electric 
and thermal transport coefhcients jointly as follows: 



(43) 



(Qx) / \aT Rt]\ -{VxT)/T 

where Zx is the electric current and Qx is the heat current, both are in the x-direction. And cr, a 
and R are the electric conductivity, the thermoelectric conductivity and the thermal conductivity, 
respectively. Following the procedures in |60|, ISlI, we can obtain these transport coefficients which 



are listed in Table HH 

From Table |IT1 we can find that both of the thermoelectric conductivity a and the thermal 
conductivity depend on the electric conductivity a and the frequency oo. Therefore, in Fig(T] and 
Figj2] we only show the numerical results for the electric conductivity a since the rest transport 
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a 


a 


H 






z 


= 1 


C2(rf+3z-5) 






P2 




m{d—l) 1 














iuT 


T 




d = 


3 


z 


_ 3 
2 


C2(<i+3z-5) 






P2 _ 
iuT 


~T 










z 


= 2 


C2(d+3z~5) 


2uj 




P2 












iuiCi i{d' 


f3z- 


5) 


iu)T 










z 
z 


= 1 

_ 3 


C2(d+3z-5) 


2a; 




P2 




7?l(d— 1) 1 




d = 


4 


iuiCi i{d- 
C2(32-1)mv'2(-- 


f3z- 
l)/3 


-5) 


iu)T 
P2 






2 








iuT 










z 


= 2 


C2(3z-1)mV2(— 






P2 












iu)C\ 






iuiT 


T 





TABLE II: The various conductivities for different z and d. 




1 2 3 4 5 10 20 30 40 




1 2 3 4 5 10 20 30 40 



FIG. 1: The real and imaginary parts of the conductivity for d = 3 with respect to various z. The left parts 
are of the low frequency regime while the right parts are of the high frequency regime. 

coefficients can be easily obtained from a. In the numerical calculations, we have scaled I = l,rh = 
1, and p2 = fJ- = I. 

Actually, in the numerical calculations, we have set the horizon very close to but not equal to 
rh, because the coefficients of the eq. (jl7p will diverge at r = r/^. From Fig{T]and FiglH we can find 
that at w = 0, the real parts of the conductivity are finite; however, the imaginary parts of the 
conductivity will diverge at cj = 0, thus from the Kramers- Kronig relations we can readily deduce 
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2 4 6 8 10 12 10 20 30 40 50 




2 4 6 8 10 12 10 20 30 40 50 



FIG. 2: The real and imaginary parts of the conductivity for d = 4 with respect to various z. The left parts 
are of the low frequency regime while the right parts are of the high frequency regime. 

that the real parts actually will develop a delta function at a; = 0. This delta function is due to 
the translational invariance of the system. These are known in previous literatures 61 1. 



For large frequencies, the expansions for a^x can be found in Table |T] in which the coefficients C2 
and C\ are functions of uj. Therefore, from Table HIl as well as Table [U we can get the approximate 
behavior of the conductivity depending on the frequency u as, 



0-d=3(w) 



/ ,0 
w , 

,2/3 



1; 



w(6 + log(w)), z = 1; 

z = 3/2; 



(44) 



a;3/2 z = 2. 



z = 

tj-/-, z = 3/2; o"d=4M ~ < 
ti;(a + log(a;)), z = 2. 

where, a and h are some constants. This large frequency behavior of the conductivities can be seen 
from the right parts of Fig. [1] and Fig. [2j 

In Figdl the real part of the conductivity will tend to a constant when uo becomes large for 
z = 1, which is similar to the previous papers 



6C 



6ll | . But the differences are in the case of z = 3/2 



and z = 2, in which the Re(cr) will depend on uj according to eq. (|44p . This is an interesting and 
new phenomenon from the viewpoint of the gauge/gravity duality, which is not observed before in 
the previous literatures as far as we know. For example, in GlJ the author argued that the electric 
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conductivity in the case of d = 3 will tend to a constant because of the dimensional analysis. 
However, here we can explicitly see that in our model for d = 3 and z > 1, a will be proportional 
to a;^(^) in the large frequency limit, where s is a function of z. This peculiar frequency dependent 
AC electric conductivity may be related to some new materials in the realistic world. Fortunately, 
in 62] the author has studied the AC conductivity for various disordered solids in (ci = 2 + 1) 
and (d = 3 + 1) dimensions both experimentally and theoretically. We found that the electric 
conductivity for d = 3 and d = 4 in our Fig(T] and Figl2] have similar behaviors to the experiments 



or the computer simulations in the large frequency limit in the paper [64]. In that paper, the author 
has proposed a kind of 'symmetric hopping model' to illustrate the large frequency behavior of the 
electric conductivities. Therefore, we expect that the Lifshitz black brane model in the present 
paper may be related to this kind of 'symmetric hopping model' from certain aspects. We will 



further report this kind of relation in another work 



In Figl2]for z = 1, the large frequency behavior of the conductivity is like u}{b + log(a;)), which 



resembles the expansions in the Appendix in 



641 ]. The arguments for the conductivity for z = 3/2 



and z = 2 are the same as those for d = 3 in Fig. [H 

IV. THE SHEAR VISCOSITY 

As we know that any interacting field theory at finite temperature in the limit of long time and 
long wavelength can be effectively described by hydrodynamics. In this section, we will compute 
the shear viscosity of the dual field theory in the low frequency limit. To do so, we need to turn on 
the transverse tensor mode fluctuation (which is in the scalar channel) of the metric Sgfj,,^ = hxy 

A. The case of iV = 1 

Let us begin with the = 1 case first, see eq.Q. Taking the mode expansion of the fluctuation 
^gxyit^f) = hxy{r)e~^^^'^^^'^ (where C, = x^-i, see the Appendix), we obtain the linearized Einstein 
equation of the xy component as, 

„ fz + d E'\ , //2+2.^2 ^2;4X 

^ + + H j ^ + (;3T2lH2 - j = 0. (45) 

which is the equation of motion of a minimally coupled massless scalar field propagating in the 
unperturbed spacetime background, where we have defined ip = hy. 

^z+d-l 

To solve eq. (j45p . it is convenient to introduce the new coordinate = p+drrr, then the boundary 
is located at u = 0, while u = 1 is the horizon. Thus the fluctuation equation taking the long 
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wavelength limit k'^ — )■ 0, becomes 



\ , X 2{z-d+l) 



u 



where = 1 — and ""' is the derivative with respect to u. At the horizon, since we are going 
to calculate the retarded Green's function of the dual field theory, we need to impose the incoming 
wave boundary condition. Thus, we set ip = [1 — u)°'^{u), then a can be determined through the 
near horizon expansion of eq. ()46p . which gives a = — p^^p^i^r^j- To obtain the solution of ^(u) in 
the full spacetime region, we can expand it in terms of cj as 

^{u) = ^oiu) + oj^iiu) + 0(^2), (47) 

and then solve the above equation order by order. Furthermore, requiring 'I'o to be regular at the 
horizon and normalizing it to be one at the boundary, as well as ^'i vanishes at the horizon, we 
find that 

IK 

then we have 

iio , , 1 + ti , 



^ = (1 - ^) 1 - , . In(^) . (49) 



To compute the shear viscosity of the boundary field theory, we need to compute the fiux factor 
J- = K^/—gg'^^(p*{u)duf{u), where K is a normalization constant related to the effective coupling 
constant of the bulk transverse graviton. Keeping to the order of 0{uj), it is straightforward to 
compute the flux factor and the retarded 2-point Green's function as 

so the shear viscosity can be obtained by the Kubo formula as 

a; " 16^G,+i;<^-i' 

then we have 

s 47r 

which satisfies the KSS bound in the Einstein gravity. 
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B. The case of > 2 

Now we consider N > 2 cases, see eq.(l8|). As we have shown in the Appendix, the equation of 
motion for {p = hy is also that of a minimahy coupled massless scalar field, which is of the same 
form as eq. ([l5]) 

^ (0 + 1^7 + —;- J (0 + [^-^j^Tfp - 747J ^(0 = 0. (53) 

Note that since /(r) has multiple zero roots and cannot be determined in general, to solve eq.(|53p. 
it is more convenient to apply the matching method in which the exact form of /(r) is not involved. 
In the near horizon region, i.e. r — r^ <^ rh, f{r) — f'{f'h){r — rh), then eq. (j53p can be simplified 

as 

'P'^r) + -^^'{u) + (^^^ - ^(r) = 0, (54) 

r-Vh \{r-rhy r-rjj 

in which 

Let's further defining f = r jr^ and taking the long wavelength limit A:^ — )• , eq. (j54p becomes 

^\f) + ^^^{f) + Y^^^^r) = 0, (56) 

which gives 

<^{f) = C'iir - \)^ ^ CA_{r - \)~^ , (57) 
in the r coordinate the solution is 

H^ir) = C'i{r - Th)'^ ^ CAir - Th)'^ , (58) 



where tt) = 2^^- The first part of (j57p or ()58p is the outgoing mode while the second part is the 
ingoing mode. To calculate the retarded Green's function, we need to adopt the ingoing mode, 
which require C3 = C3 = in eq.(|57p and eq. (|58p . In the low frequency limit, eq.(|58p can be 
expanded as 



^(0 = C4 1^1 - ^ ln(r - r,) + 0(u;^)J , (59) 
In the near region, r/jO; < ra; ^ 1, then in the A;^ — t- limit, eq. (|53p reduces to 

^"(r) + + ^) ^'(r) = 0, (60) 
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which can be solved as 



= / + (61) 



Note that in the near horizon hmit r — )■ r/i, eq. (|6ip can be simphfied as 



f^'^) - I T-dTz'^'^ + C6 = ln(r - Th) + C6. (62) 



^ -L (^^1 = ^ 

/'(r/.)(r-r,)rf^ /'(r,,)< 
While in the large radius limit, /(r) — ?• 1, then eq. (l6ip becomes 



^(0 - / -^Ar + C6 = + C6, (63) 

In the outer region <^ I <^ r, f'{r) — t- 0, /(r) — ?• 1 and again we taking — )■ 0, then eq. ()53p 
becomes 

/(O + ¥^'(0 + ^;To^V(r) = 0, (64) 



r 



in the n = 1/r coordinate, eq. (j6l]) can be changed to 



ant its solution is 



(p[u) =U 2 ( C7J_A+ ( I +C8Ja+ ( ) ) , (66) 

where 

cj = cj{2z) 22 {V+''uj) 2z r and cg = C8(22;) (I'-^'^uj) 2z T 



2z J ^ ' ' ' \ 2z 

in which, C7 and cs are certain constants while A+ = d + z — 1 is the conformal dimension of the 
operator dual to the massless scalar field in the bulk. Again, in the low frequency limit, eq. (|66p 
can be expanded as 

(^(r) = cr (1 + 0(0.2)) ^ ^-/^^ 0^ ^ u;i+^r-^+ (l + ©(a;^)) ^ (g^) 

The condition for matching the solutions in these three regions is rh < r <^ oj^^ . Comparing 
the eq. (l62]) with eq. (|59]l we get that 

C4 = ce and — iojC4, = '^^ d-i ' (^^) 



While the matching of eq. ()63]) with eq. (l67P gives 
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Namely, the coefficients in tfiese tliree regions are related by the following relations 

(l+z)A+ /2\^ i,d-l P+V^"^ ^ . 

cy = ce = C4 and 03/ ^ - uj ^ = C4- iu;. (70) 

\z J d + z — 1 

Furthermore, the normalization condition requires that (p{r) is normalized to be 1, namely, C4 = 1. 

Consequently, the asymptotic solution at low frequency limit becomes 

y,(r) = (r-r,)— il + ^^^-—^iu; + 0{u^)\ . (71) 

After eliminating the divergent terms, the dominant part of the radial flux of the scalar field at 
the boundary is 



J- = K^g"-^*{r)dMr)\r 
1 

-uj + 0{lo-'), (72) 



where K = 1/{32ttGii^i) is the effective coupling constant of the scalar field (p{r), then the retarded 
Green's function is 

Gnik) = -2T{k,r)\r-,oo, (73) 

and the shear viscosity is calculated from the Kubo formula 

IuiGr 

r] = — hm 

h (74) 



Therefore, the ratio of the shear viscosity to the entropy is 

; - h 

which gives the same value as that of the Lifshitz black brane with only one gauge field. 

The result indicates that the additional background U{X) gauge fields do not alter the KSS bound 
of the boundary fluid, though they do contribute to the shear viscosity and the entropy density, 
respectively. 



V. CONCLUSIONS AND DISCUSSIONS 



In this paper, we studied the model of strongly coupled non-relativistic quantum field the- 
ory with multiple C/(l) gauge fields near the Lifshitz fixed points, in the framework of the non- 
relativistic gauge/gravity duality. By considering the linearized perturbations of bulk gravitational 
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and gauge fields, we solved the equation of motions for gauge fields with backreactions (shear 
channel) and the bulk transverse graviton (scalar channel). For the N = 2 case, we derived the 
renormalized second order effective action and systematically calculated the electric, thermal and 
thermoelectric conductivities of the dual non-relativistic quantum field theories with respect to 
various d and z. Specifically, we found the novel frequency dependent power law behavior of the 
AC electric conductivity in the large frequency limit when d = 3 and z > 1. From the knowledge of 
the condensed matter physics, we expect that our model provides a holographic description of the 
'symmetric hopping model' in some sense. The argument goes to the case of d = 4 as well, we will 
report the further relationship between the Lifshitz black brane and the hopping conductivities in 
another paper elsewhere. In addition, when taking the limit of long wavelength and low frequency 
in the generic N cases, we also showed that the ratio of shear viscosity to entropy density of the 
dual boundary fiuids still satisfies the KSS bound derived in the Einstein gravity. 
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The Einstein-Maxwell-dilaton theory with multiple U{1) gauge fields that we are considering 
has the action 



Acknowledgement 



Appendix A: Linearized Perturbations of the Gravitational Theory 



1. Einstein-Maxwell-dilaton theory 



I = 



1 




a=l 



) 



(Al) 



its Einstein equation is 




2A 



.A 



1 



9(11/ 



(A2) 



d-1 



au 



2{d-l) 
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Let us make the metric ansatz to be a d + 1 dimensional black brane solution as 

ds'^ = Hi{r) {-f{r)dt^ + dx'dxi) + H2{r)dr^, (A3) 

where its outter horizon is located at f{rh) = 0. 

Consider the small metric fluctuation caused by some external perturbation 

jS'J ^9f.u = 9i!'J+6g,., (A4) 



the Christoffel symbol is 

pA _ p{0)A _L r-pA 
„Aa 

= rg]^ + ^— {V^dg^, + V,%, - Vc^Sg^,) , (A5) 

when taking the linear order perturbation of the metric, i.e. dg^^i, = /i^^y, the Christoffel symbol 
can be expanded up to second order of h as 

= rw^ + r«^ + r(i)\ (A6) 



where 



{0)Aa 

rW>^ — 1 (V h -I-V7) — V/) 1 

rJS^ = — ^ (V;,/i«^ + V - v« V) ■ (A7) 



Note that under the first order variation, the Ricci tensor varies as 



Ri^)+RW+Ri^l (A8) 



where 



^ ■p(l)a _-p(l)a 



^(V"V^/i«^ + V^V^K^) - - ^V^V^/i, (A9) 



and 



7,(2) ^ U2)a _ U2)a 



■^TT- (^mV + Vz./1;3m - V/3 V) + ^^V^/i„^, (AlO) 
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the first order and second order Ricci scalars are 

d-1 



2A 

V"V^/i„^ -nh- - — -h, (All) 



and 



^(2) ^ ^{0)A..^{2) _ J^^.u-^{l) 



-V«/i"^V^ V + ^^Va/. + ^^V^ V> (A12) 



respectively. 

Then the linearized Einstein equation is 

N 



a=l ^ ^ ' ' 

When there is only transverse gravitational fluctuation hxy = hxy{r)e~^^^^'^^^ , where C, = x^-i is 

(0)x 

the d — 1-th spatial coordinate. Using the Rx component of the zeroth order eom, i.e. 

a=l \ \ / / 

then eq.dSS]) becomes 

_ ^Uip + ( - -DhA ^ = -^H.ip - - y e^"'^ , ^ FlH^^ = 4% (A14) 
2 ^ \2HiH2 2 d-1 ^ 2^ 2(d-l) " ^ xxH^ \ I 



which gives the eom of minimally coupled massless scalar field (p = hy 

Hi 1 
2 



^in^ = -^^1=9^ [^/^^g^'^'-'d,^^ = 0. (A15) 



2. Gauge field perturbation with backreaction 

To compute the conductivities of the dual field theory, we need to turn on the gauge field 
perturbation along the spatial direction, this gauge field perturbation will in turn induce the hu 
off-diagonal part of the background metric perturbation since hu and both belong to the vector 
mode fluctuations. Without loss of generality, we choose SAa^ = 6xaax{r)e~^''^^~^'^'"^ , which induces 
the corresponding metric perturbation is Sg^i, = htx{r)e~^'^^~^^^'^ . Then the linearized Einstein and 
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Maxwell equations are obtained by making the combined diffeomorphism and gauge variations to 
the original equations, namely 



+ ^E^^""^ (^-A^^i - ^^(^fla,^ , (A16) 



where 5e means the diffeomorphism transformation while 5^ indicates the gauge field transformation 
that obeying the following relations 

SeQiiu = CtQiiu and (^x^^M = (A17) 

In the linear order perturbation, the nonvanishing components of the first order Ricci tensor are 
Rxt = Rtx and Rxr = Rrx- Then the linearized Einstein equation are 

^^'^ - J^l^^' = \ E [a'^'^'-Qra^^QrA^t - ^ (A18) 

together with the xx component of the zeroth order Einstein equation, eq. (|A18p becomes 
1 



1 ^ 

Y.e^^'^a'^^A',,, (A19) 



a=l 



and 



N . N 

^rt^ = E e'^'a^'^^'d^a^^Kt = f E e'^'^g^'^^'aaxKu (A20) 

a=l a=l 



which gives 

h'tx - 'irhtx + Yl ^^^^^axA'at = 0, (A21) 

where ""' indicates d^.. 

In addition, the linearized Maxwell equation is obtained by 



(y^e^'"^9^°)^°5^°^'^^(5aaa;3 - dpaaa)^ = 0, (A22) 
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its nonvanishing components are 

+ dr [y^ge^^^g^'^^'^g^'^^^^draa.) = 0. 
In the black brane background eq. ()A3p . the above equations become 

When taking the ansatz Hi = 6^, H2 = 1/^, /(r) = ^e~^/6^ and 70(0) = e'^"'^ in eq.(j4]), 
eq. ()A2ip and eq. (|A23P change into 

26' ^ 

- -yhtcc + lai(p)Ktaax = 0, (A24) 

a=l 

and 

In the hnear order perturbation of the metric and the gauge fields, the bulk action can also be 
expanded into second order as 

/ = /(o) +/(!)+ /(2), (A26) 

where the zeroth order action is 

and the first order action is 

= Y^;^ d^x^^n^ l^'^K - ^''^ - E , (A28) 

which are purely surface terms when the bulk eoms are satisfied (on-shell condition), where is 
the unit normal vector of the hypersurface S. 
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While the second order action is 

_£ Vp^""?^ f f(i) 7?(i)a''^ _ 47?(i) I 77(0) 17(0) 

a=l 

where jC^*^^ and C^^^ are respectively the first and second order Lagrangian densities in /'^'^' and 
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